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Convection in binary fluid mixtures: A model of four 

coupled amplitudes 
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We consider the bifurcation scenario that is found in Rayleigh-Benard convection of binary fluid mixtures 
like ethanol-water at positive separation ratios and small Lewis numbers leading to a bifurcation sequence 
of square, oscillatory crossroU, stationary crossroU, and roll patterns. We propose a system of four coupled 
amplitudes that is capable to model all important properties. 
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1 Introduction 

The Rayleigh-Benard system ^ |2] is one of the 
most popular model systems to study spontaneous 
structure formation emerging via instabilities. The 
system consists of a fluid layer confined between 
two parallel plates perpendicular to the direction of 
gravity, and exposed to a vertical temperature gra- 
dient that is established by keeping the upper and 
lower plate at two different temperatures Tq — AT/2 
and To -|- AT/2 respectively. 

Rayleigh-Benard convection in one-component 
fluids at small and moderate temperature differ- 
ences is well understood, mainly thanks to the 
work of Busse and co-workers in the 70s and 80s 
|S1 in El El • The bifurcation behavior becomes more 
complex in binary mixtures. Here, the concentra- 
tion enters as an additional relevant dynamic field. 
Concentration differences are created in the pres- 
ence of a non-vanishing Soret effect, a nonequilib- 
rium effect that describes the driving of concentra- 
tion currents by temperature differences. 

The concentration couples back into the system's 
dynamics via variations that determine the buoy- 
ancy force density. When the lighter component of 
the mixture is driven into the direction of higher 
temperature, thereby enhancing the density gradi- 
ent and further destabilizing the layer one speaks of 



a positive Soret effect, and of a negative Soret effect 
in the opposite case. 

Convection in binary mixtures shows a very rich 
bifurcation behavior. For negative separation ra- 
tios convection rolls bifurcate backwards and time- 
dependent patterns, namely traveling waves and 
standing waves appear. The bistability of convec- 
tive patterns and the ground state leads to the ex- 
istence of localized structures and fronts. 

In the case of positive separation ratios on which 
we will focus in this paper the two-dimensional roll 
structures are often replaced by squares at onset, 
but become the preferred pattern at higher temper- 
ature differences. For intermediate AT two types of 
crossroU structures can exist, a stationary and an 
oscillatory type that transfer stability from squares 
to rolls. 

There are few-mode models that satisfactorily 
explain the bifurcation properties for negative Soret 
effects [71. However, for positive Soret effects a 
model with only a few degrees of freedom, that is 
able to properly describe rolls, squares and cross- 
rolls seems to be lacking. 

In this paper we will present a system of four 
equations that describe all important properties of 
the square - crossroU - roll transition for the case 
of a positive Soret effect. 

The paper is organized as follows: after this in- 
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troduction we discuss in Sec. 2 the Rayleigh-Benard 
system and the bifurcation scenario for positive 
Soret effects in more detail. In Sec. 3, the main part, 
we present and discuss a model system that shows 
a similar behavior but incorporates an additional, 
artificial symmetry that makes it easier to identify 
the stationary solutions but leads to some qualita- 
tive disagreements with the full system. In Sec. 4 
we add some terms that break this symmetry and 
demonstrate that then these differences disappear. 
We summarize our findings in Sec. 5 and provide an 
outlook. 

2 Bifurcation scenario 

The basic equations that couple convection velocity 
u, temperature 6, concentration c and pressure P 
read in a dimensionless form [S] 



{dt + u-V)u 

{dt + u-V)c 
Vu 



-VP + c 
Ru-, + V 



{6 + c)e, + V^u 







The fields are given as deviations from the qui- 
escent conductive ground state. The direction of 
gravity is the negative z direction. See ref. |S] 
for details. There appear four different parameters 
in these equations: the Rayleigh number R is the 
control parameter, the dimensionless temperature 
difference between upper and lower plate. In pure 
fluids, the critical Rayleigh number of the onset of 
convection is -Rc(V' = 0) = 1707.76. One defines a 
reduced Rayleigh number r = R/Rc{^ = 0). The 
Prandtl number a is the ratio of the timescales of 
heat and momentum diffusion. The Lewis number 
L on the other hand is the ratio of the heat and 
concentration timescales. The separation ratio fi- 
nally measures the strength of the Soret effect. In 
the conductive ground state it determines the ratio 
of the density gradients due to concentration and 
temperature. 

In typical liquids it is o" ~ 10. Concentration 
diffusion is slow in ordinary liquid mixtures; a typ- 
ical value for L, say, for ethanol-water is 0.01. ip 



varies with the mean concentration and mean tem- 
perature and can be both positive and negative in 
ethanol- water. 

For positive ■0 three different bifurcation scenar- 
ios exist at small and moderate r. (i) For very small 
-0 and relatively large L the Soret generated con- 
centration variations are small and easily diffused 
away. The mixture then still behaves qualitatively 
like a pure fluid. The first convection structure to 
be observed above onset takes the form of parallel, 
stationary rolls with alternating direction of rota- 
tion. Near onset, Uz, 0, and c can be described by a 
function Acos{kx)f{z) with k being the wavenum- 
ber of the rolls and f{z) being different for the three 
fields and defining the critical vertical profile. 

(ii) When ip is larger and/or L is smaller the mix- 
ture begins to behave qualitatively different. The 
rolls are replaced by a square pattern near on- 
set that can be written linear combination 
of equally strong rolls in x- and y-direction, i. e., 
[Acos{kx) + Bcos{ky)] f{z) with A = B. At higher 
r however the advective mixing is strong enough to 
equilibrate the concentration field outside of some 
small boundary layers at the plates and between the 
rolls. The mixture again behaves like a pure fluid 
and squares lose their stability to rolls. The stabil- 
ity is transferred via a crossroll branch, a structure 
that, like squares, can be understood as a superposi- 
tion of two perpendicular sets of convection rolls. In 
the case of crossrolls, however, the two amplitudes 
A and B are not equal and the crossroll structure 
therefore lacks the square x <-> y symmetry. Like 
squares and rolls, these structures are stationary. 

(in) At even smaller L, but easily realizable with 
liquid mixtures the scenario is again different as 
shown in Fig. ^ There the squared amplitude of 
the leading velocity mode of the x-roll set is plot- 
ted versus the reduced Rayleigh number r. Solid 
(open) symbols denote stable (unstable) structures. 
Squares (square symbols) are again stable at small 
r but lose their stability against another type of 
crossrolls, an oscillatory one (lines without sym- 
bols). The bifurcation point where the stationary 
crossrolls (triangles) emerge still exists but at higher 
r, where the oscillatory instability has already taken 
place. In the oscillatory crossroll branch the ampli- 
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FIG. 1. Bifurcation scenario of the fuh hydrodynamic 
field equations for a = 27, L = 0.0045, ip = 0.23. The 
wave number of the structures is fc = tt. 

tude of the x- and y-roll set varies periodically in 
counterphase to each other such that when at one 
time the x~rolls are stronger, half a period later the 
y-rolls will be the stronger set. The two branches 
in Fig. n show the maximal and minimal ampli- 
tudes of these time dependent structures. The two 
stationary crossroll branches represent the two dif- 
ferent types: In the upper branch the x-roll com- 
ponent is the dominant one and their amplitude 
grows whereas the amplitude of the y-roll compo- 
nent gets smaller. When the latter component van- 
ishes the crossroll ends at a bifurcation point on the 
x-roll branch (circles). The lower branch belongs 
to the type of crossrolls where the y-roll compo- 
nent is the dominant one. Now the plotted x-roll 
component gets smaller until the structure ends on 
the y-roll branch, i. e., on the axis in Fig. ^ since 
there only the strength of the x-rolls is plotted. The 
stationary crossrolls gain stability when the oscilla- 
tory crossrolls vanish and transfer the stability to 
the rolls. The existence of all these structures have 
been demonstrated in experiments [HI IIUI ITT] . 

The time-dependence of A and B is plotted for 
two oscillatory crossrolls in Fig. [3 It shows the 
main characteristics of these structures near the be- 
ginning and the end of their r-range of existence. 
While the oscillation is very harmonic at smaller r 
it becomes very anharmonic at higher r. The fre- 
quency decreases. Finally, the oscillatory crossrolls 
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FIG. 2. Oscillatory crossrolls, squared amplitude of the 
X- and y-roU component for the parameters of Fig. ^ 
and r = 1.08 (top), r = 1.13 (bottom). 

end in an entrainment process in a very narrow i — 
interval jl2j . 

The crossroll structures exist around r = 1, i. e. 
the critical Rayleigh number of the pure fluid. In the 
same region the branches of rolls and squares bend 
upwards. Before and after this transition interval 
the curves are almost linear but with very differ- 
ent slopes. These two regimes at small and large r 
are called the Soret and the Rayleigh regime respec- 
tively. In the Soret region the convection can be well 
explained by an amplitude equation model includ- 
ing only cubic terms. In the Rayleigh region on the 
other hand the concentration field has been equi- 
librated and the convection amplitudes agree well 
with those of the pure fluid. In the Rayleigh region 
slightly above r = 1 the pure fluid still obeys its 
own amplitude equation and so does therefore in- 
directly also the mixture ^H]- We will make use of 
the fact that the bifurcation behavior of the mixture 
can in fact be described by two different amplitude 
equations in two different r-intervals in the next 
section. 

3 The simplified model 

In this section we want to present a flrst attempt to 
explain the bifurcation scenario as described above 
by means of a simple model. The full picture cannot 
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be captured by simple amplitude equations derived 
in an expansion around the onset in combination 
with a multi-scale analysis from the basic equations. 
Although a sequence of stable squares, forward bi- 
furcating stable (stationary) crossrolls, and finally 
stable rolls may be described with a system of two 
coupled quintic amplitude equations [H], the oscil- 
latory crossrolls are not contained in such a model. 
When A is the real amplitude of x-rolls, referring to 
a critical mode A cos{kx) f (z) , and B is the corre- 
sponding mode in y-direction then a system of two 
coupled amplitude equations must read for reasons 
of symmetry 

dtA = fiA- F (^A^,B'^^ A , 

dtB = fiB - F {^B'^ , A^^ B , (2) 

with some arbitrary function F. F is real since the 
critical perturbations lead to stationary patterns. 
For squares, it is ^ = and the Jacobian of the 
right hand sides of the above equations will be sym- 
metric. Therefore, complex eigenvalues that would 
mark the appearance of oscillatory crossrolls are not 
possible along this branch. 

In order to reflect the existence of the Soret and 
Rayleigh regime we construct a model of two cou- 
pled sets of amplitude equations. We define two 
different amplitudes Ai and A2, both describing 
rolls in x-direction of different nature: solutal rolls 
that exist in the Soret region due to the desta- 
bilizing effect of the initial concentration gradient 
and which are described by a mode Ai cos{kx)f{z) 
and the thermal rolls A2Cos{kx)g{z) that exist in 
the Rayleigh regime environment of almost homo- 
geneous concentration and which are driven by the 
thermal stress alone. Bi and B2 are the correspond- 
ing amplitudes for the rolls in y-direction. 

The most simple model would consist of two sets 
of equations as in ((2)), including only cubic terms 
to couple the A- and i3-rolls and additional cubic 
terms to couple the index 1 and index 2 structures: 

dtAi = fiiAi 

-AiihAl + ciBf + diAl + eiBi) , 

dtA2 = /U2^2 



-A2{h2Al + C2Bl + d2Al + e2Bl) ,(3) 
dtBi = niBi 

-BiihBf + ciAl + diBi + eiAl) , 

dtB2 = /U2-B2 

-B2{b2Bi + C2Al + d2Bf + e2Al) . 

These equations do not define an amplitude equa- 
tion model in a strict sense since they are not meant 
to result from a controlled expansion of the dynam- 
ics near the bifurcation threshold of the full system. 
Instead, in order to describe the actual convection 
behavior we have two bifurcation points here, as it 
is displayed by the two different control parame- 
ters /ii and fj,2- These parameters are not indepen- 
dent however. Since the only control parameter in 
the Rayleigh-Benard system is r, the parameters fii 
must be functions of each other. We set 

fJ-i= fJ- , iJ,2 = afi - fio . (4) 

The two bifurcation points are at /xi = and at 
fi2 = 0. Since for L <^ 1 solutal rolls grow on a 
much slower timescale than thermal rolls one has to 
have a ^ 1. When = — 1 refers to r = it has 
to be /X2 = at ;U » 1 since the onset of convection 
in mixtures lies at much smaller r than the onset 
r = 1 in a pure fluid. For the plots we set a = 10 
and /io = 100. That means that the timescale of 
the solutal rolls is assumed 10 times larger than the 
one of the thermal rolls. The latter appear only at 
H > 10, or 11 times the critical Rayleigh number of 
the solutal rolls. The ratios found in the numerical 
investigation of the full system are typically even 
larger for realistic parameter combinations. 

We assume all parameters 61,62, •••62 to be posi- 
tive. A different sign of any of these would lead to 
unwanted features in the bifurcation scenario. Since 
a ^ 1 we can assume that for any two parameters p, 
q € {61, 62, ...62} the inequality p < aq holds. This 
assumption will simplify our discussion below. We 
can set 61 = 62 = 1 by normalizing the amplitudes 
accordingly. 

The equations ^ show several symmetry prop- 
erties. They are invariant under an exchange of A 
and B amplitudes, i.e., under exchange of x and 
y. The absence of quadratic terms guarantees the 
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preservation of the upflow-downflow symmetry of 
the convective fields. A change of the sign of only 
Bi and B2 or only Ai and A2 corresponds to a shift 
of the X- or y-rolls by half a wavelength. 

Furthermore, it is even possible to change the sign 
of only one amplitude. However, this symmetry is 
artificial. The real convection patterns must be con- 
sidered nonlinear superposition of the solutal 
and thermal patterns. Changing, e. g., in x-rolls 
the direction of rotation of the solutal contribution 
while leaving the thermal one unchanged is certainly 
not a symmetry operation. Nevertheless we will add 
terms that break this symmetry of eqs. ^ only in 
the next section. 

Although the bifurcation scenario of the simpli- 
fied model © disagrees qualitatively with that of 
the full hydrodynamic field equations in some points 
it nevertheless facilitates to understand some other 
features. Furthermore, the simplified model easily 
allows to get analytical solutions for the stationary 
structures: Every single mode can be set indepen- 
dently to zero. The remaining equations can in the 
stationary case then be rewritten as a linear system 
with the squares of the nonzero amplitudes as vari- 
ables. The unique solution is a linear function of 
;U, and physically relevant in a /^-interval where all 
amplitude squares are > 0. We can conclude that 
there are at maximum 16 different stationary solu- 
tions of eqs. © because there are 2^ possible ways 
to set some amplitudes to zero. 

3.1 Thermal and solutal rolls 

The first kind of structures we discuss are x-rolls 
with Bi = B2 = as solutions of the subsystem 



dtAi 



HiAi - Ai [biAj + diAl 



dtA2 = fi2A2 - A2 (b2Al + d2Af^ 
Since 61, 62 > 0, both pure solutal or ^i-rolls 



(5) 



Ml 



Ao 



(^i = /i>0) , (6) 



and pure thermal or -rolls 

Al = ^ , Ai = (^2 = a/i - /xo > 0) , (7) 
02 



exist as forward bifurcating solutions. They are 
plotted in Fig. |31for 61 = 62 = 1 and a, /io as above. 



100 



< 




FIG. 3. Bifurcation of pure solutal {Ai) and pure ther- 
mal {A2) rolls of the model eqs. ©■ Parameters are 
a — 10, ^0 ~ 100, and 61 = &2 = 1- 

As the only convective solution for small fi, the 
solutal rolls are stable at onset. But the linear sta- 
bility analysis shows that they lose stability in fa- 
vor of the stronger thermal rolls at larger /i. The 
stability is transfered via a third kind of rolls, the 
intermediate or ^i^2~i'olls given by 

62M1 - diH2 



^2 



6162 
bifJ-2 



did2 
d2fJ-i 



bib2-did2 ■ 

These structures exist if there is an /^-interval where 
both right-hand sides are > 0, what is always the 
case. If 6162 — did2 > they bifurcate forward from 
solutal to thermal rolls. 

A bifurcation diagram showing all three kinds of 
solutions and also their stability properties within 
the restricted system © is shown in Fig. 21 Plot- 
ted is only the interesting |U-range where stability is 
transferred. The stable branches already look like a 
caricature of the roll branch in the full system with 
different slopes in the Soret and Rayleigh regions. 
However, there are also important qualitative dif- 
ferences. In the full system, the roll branch consists 
of one continuous structure, without the two bifur- 
cation points that are found in the model system. 
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These bifurcation points will disappear when we 
break the artificial symmetry present in the model 
system. 
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FIG. 4. The three types of stationary rolls solutions 
of the model eqs. (jSjl: Solutal rolls (Ai), thermal rolls 
{A2) and intermediate rolls {A1A2) plotted as Af + 
versus [i. Sohd (dashed) lines denote stable (unstable) 
structures. Parameters are d\ = di= 0.88 and otherwise 
as above. 



3.2 Thermal and solutal squares 

A solution has square symmetry when A\ = B\ 
and A2 = Bi- Analogous to rolls there are three 
different types of squares, namely solutal or A\B\- 
squares 



B 



hi+ci 
thermal or ^2-B2-squares 



^2 = 



^2 



Bi 



A, =B, 



, 



62 + C2 

and intermediate AiBiA2B2-squaTes 

(62 + C2)Ati - {di + ei)ii2 



(9) 



(10) 



Aj = Bf 



Al = Bl 



A 

(61 + Ci)/i2 - (^2 + e2)/il 



,(11) 



where 

A = (6i+ci)(62 + C2)-(di + ei)((i2 + e2) . (12) 



Under the general requirements stated above the in- 
termediate squares exist and they bifurcate forward 
from solutal to thermal squares if A > 0. The bi- 
furcation branches of the different square and roll 
structures are plotted in Fig. [3 
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FIG. 5. Square and roll structures of the model eqs. 
Squares (rolls) are represented by thick (thin) lines. Solid 
(dashed) lines denote stable (unstable) structures. Bi- 
furcation points of further interest are marked by circles. 
Parameters are c\ — 0.96, C2 — 1.04, ei = 0.68, 62 = 1.28 
and otherwise as above. 

The stability properties of the different types of 
squares against each other are analogous to those 
of the rolls. But these are not all types of possible 
instabilities. Since our model has four modes four 
eigenvalues have to be investigated for all structures 
found so far. We will briefiy discuss the respective 
results. 

Some instabilities are connected to the bifur- 
cation points that are known already. For the 
squares only those instabilities that can break the 
{x y)-symmetry are still of interest. It turns 
out that solutal squares are stable against (x ^ 
y)-antisymmetric ((^74i, (^i?i)-perturbations every- 
where if and only if c\ < bi. This inequality is 
satisfied in the example case and therefore squares 
are stable at onset as desired. Similarly, thermal 
squares are stable against a type of {6A2, 6B2)- 
perturbations only if C2 < 62. This is not the case 
and consequently the thermal squares have been 
marked as unstable in Fig. [5J 

These requirements have also consequences for 
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the stabihty of rohs: Solutal rolls, although sta- 
ble against other x-roll structures for small fi are 
unstable against (5i?i-perturbations if and only if 
ci > bi as required for solutal squares to be sta- 
ble. On the other hand, thermal rolls are stable 
against (5i?2"Perturbation when C2 < &2 and thermal 
squares are unstable. Therefore, only the thermal 
rolls are marked as stable in Fig. [51 

There are two other types of instability of rolls 
that have to be discussed, namely the instabil- 
ity of solutal ^i-rolls against (5i?2~perturbations 
and vice versa, of thermal A2~rolls against 5Bi- 
perturbations. Both instabilities generate new bi- 
furcation points that have been marked by large cir- 
cles in the figure. Their position depends on the di 
and Cj. Due to the fact that 62 > d2, the new bifur- 
cation point on the solutal roll branch lies at higher 
^ than the bifurcation point to the intermediate 
rolls, where the solutal rolls already lose stability. 
Since ei < di, on the other hand, the new bifur- 
cation point for thermal rolls lies on the otherwise 
stable part of the thermal roll branch, destabilizing 
it towards the region of smaller ^ already before the 
bifurcation point towards the intermediate rolls. 

Placing the new bifurcation points this way no 
new instabilities are to be expected for the in- 
termediate rolls. They inherit stability against 
(^i?2"Perturbations and instability against 5Bi- 
perturbations from the solutal rolls and transfer 
these properties to the thermal rolls. 

It should be pointed out that by choosing the pa- 
rameters accordingly it is also possible to stabilize 
the rolls in the whole /i-interval against perturba- 
tions in y-direction, which is the scenario that can 
be found at large L and small ■0 where the square 
patterns are always unstable. 

Concerning the square structures we have to re- 
mark that at the bifurcation point where the inter- 
mediate squares branch off from the solutal squares 
in fact not one but two eigenvalues go through zero 
for symmetry reasons, namely both 6A2~ and SB2- 
perturbations. Two types of crossroll branches ap- 
pear here plus the intermediate squares as nonlin- 
ear combination of the two. The same happens at 
the bifurcation point of the intermediate and ther- 
mal squares. Here also two new patterns emerge, 



connected to growing 5Ai- and (5Bi-perturbations. 
These bifurcation points are also marked by circles 
as a reminder of the new branches that still have to 
be discussed. 

3.3 Crossroll structures 

We have already pointed out above that at maxi- 
mum 16 different stationary structures exist within 
the model ©. The most simple one is the ground 
state where all amplitudes vanish. Furthermore we 
have already discussed the Ai-, A2-, and A1A2- 
rolls that have their counterpart in y-direction, i. e., 
the Bi-, B2-, and i?ii?2 "rolls. Together with the 
A1B1-, A2B2~, and ^i-Bi^2-B2-squares these are 
10 structures so far. The remaining six are crossroll 
structures, i. e., three-dimensional patterns lacking 
the (x <-> y)-symmetry of squares. 

The first two of these structures are the A1B2- 
crossrolls and the Si^2~crossrolls that can be 
mapped onto each other via the (x y)-symmetry 
operation. The Aii32-crossrolls connect the Ai~ 
rolls and the i?2-i'olls at the newfound bifurcation 
points for the roll structures. Note that i?-roll 
branches would lie on the //-axis in the plots since 
only Ai + A2 is plotted. The stability analysis re- 
veals two more bifurcation points on these branches 
as it has to be since they connect the stable thermal 
rolls to the twofold unstable solutal rolls. 

The next pair are the A1B1A2- and A1B1B2- 
crossrolls that branch off at the connection between 
solutal and intermediate squares and end up at one 
of the bifurcation point of the B1A2- and A1B2- 
crossrolls, respectively. 

The A1A2B2- and i?i^2-B2-crossrolls finally 
emerge at the intermediate squares /thermal squares 
bifurcation point. These too end up at the A1B2- 
and i3iA2-crossroll branches, at the other two re- 
maining bifurcation points. 

That completes the description of the station- 
ary structures. The results are shown in Fig. |21 
together with the result of the full stability analy- 
sis. The unimportant part of the solutal and ther- 
mal roll and square branches have been cut off for 
better visibility. At small /i solutal squares are sta- 
ble. They lose their stability where two types of 
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crossroUs {A1B1A2, A1B1B2) and the intermediate 
squares branch off. The intermediate squares gain 
stabihty whereas the crossrolls are unstable. Tak- 
ing only stationary structures into account the in- 
termediate squares remain stable until a second bi- 
furcation point, where again two crossroll branches 
emerge {A1A2B2, B1A2B2) that are now stable. 
The thermal squares on the other hand do not gain 
stability. Next, the stable crossrolls meet a third 
pair {A1B2, -61^2) that finally transfers stability to 
the thermal rolls. The bifurcation diagram looks 
qualitatively different for other parameter combina- 
tions but the inequalities that have to be fulfilled 
to produce a diagram as in Fig. IHl can be derived 
easily. 

The picture is still not complete however. Al- 
though all stationary instabilities have been iden- 
tified now, there might be more bifurcation points 
where time dependent patterns emerge. From the 
otherwise stable branches only the mixed squares 
and the BiA2B2~ and ^iA2-S2-crossrolls can ex- 
hibit oscillatory instabilities. For the given parame- 
ters such instabilities do indeed occur. They desta- 
bilize the parts of the branches between the circles 
in Fig. ini One can easily guess that this is the re- 
gion where oscillatory crossrolls can be found. But 
we want to discuss these structures in more detail 
only in the full model. 
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FIG. 6. Stationary crossroll solutions of eqs. Thick 
(thin) lines represent crossrolls (rolls and squares). Solid 
(dashed) lines denote stable (unstable) structures. Os- 
cillatory instabilities occur at the circles. 



4 The full model 

The bifurcation diagram we found does already 
show some of the properties of the scenario we 
want to explain but turned out to have unwanted 
bifurcation points. We will now show how in- 
cluding terms that break the artificial extra sym- 
metry leads to a much better model. There are 
four more cubic terms preserving the symmetries 
(^1,^2) ^ {Bi,B2), {Ai,A2) ^ (-^1,-^2), and 
{Bi,B2) ^ {-Bi,-B2). When denotes the 
old right hand sides of Q we enlarge the model as 
follows 



dtAi 



dtA2 



dtBi 



dtB2 



-A2{f3iAl + 71^2 + ^1^12 + eiBl) , 



-Ai{(52Al + 72^1 + 52AI + e2Bl) (13) 



-B2{(5iBl + 71^2 ^ ^^^2 ^ ^^^2) ^ 



-Bi{l52Bl + 72^1 + 52BI + e2Al) 



The new equations do not allow a sign reversal of 
individual amplitudes anymore. However, the pres- 
ence of the /3j, 7i, and alone still allows pure solu- 
tal or thermal stationary roll solutions. This leads 
again to unwanted bifurcation points. On the other 
hand, when (5j 7^ such solutions become impos- 
sible. We will therefore discuss only these two and 
set the other six to zero. Thus the model we discuss 
here is 

dtAi = niAi - 5iAl 

-Ai{biAj + ciBf + diAl + eiB^) , 
dtA2 = 112A2 - 52AI 

-A2{h2Al + C2Bl + d2A\ + e2Bl) (14) 
dtBi = fiiBi-5iBl 

-BiihBf + ciAl + diBl + eiAl) , 

dtB2 = f^2B2-S2Bf 

-B2{h2Bl + C2Al + d2Bl + e2A\) . 

It is instructive to first take a look at the im- 
pact that the two new symmetry-breaking terms 
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have on the eigenvalues of the important branches 
as they emerge from the linear stability analysis. 
Fig. [3 shows the real parts of the three most impor- 
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FIG. 7. Eigenvalues (EVs) along the important parts of 
the different roll and square branches for 5i = 62 = 0.01. 
Dots denote the simplified system with 6i =Q. 

tant eigenvalues for both rolls and squares for small 
values 5i = S2 = 0.01 and 6i = in comparison. 
Plotted are the intermediate structures in their ^- 
interval of existence and the solutal and thermal 
structures below and above, respectively. In the 
upper plot, one can see how one eigenvalue of the 
solutal and thermal rolls becomes zero in the sim- 
plified system at the points where the intermediate 
rolls emerge. The referring curves would cross the 
zero axis here if the picture would not switch to the 
intermediate rolls. For 6i ^ these bifurcations be- 
come imperfect. The curves become smooth and do 
not touch the zero axis anymore such that the two 
bifurcation points vanish. The respective parts of 
the solutal, intermediate, and thermal roll branches 
fuse into one single branch. On the other hand, 
the bifurcation point to the crossrolls at /z ~ 11.7 



where the rolls finally become stable does not go 
away when <5j 7^ 0. 

In the lower plot one can observe a similar be- 
havior for the square eigenvalues. In the simplified 
system the solutal rolls are stable until one eigen- 
value crosses the zero axis and the picture switches 
to the now emerging intermediate squares. This bi- 
furcation and also that between intermediate and 
thermal squares becomes again imperfect for 5i ^ 0. 

For 6i = 0, at the beginning of the intermediate 
square interval two eigenvalues approach each other 
on a very small ^-interval such that the curves look 
discontinuous here. When they meet they form a 
complex pair and only the real part is shown now. 
This complex pair becomes critical in the middle of 
the interval and gives rise to the emergence of the os- 
cillatory crossrolls. After the pair separates rapidly 
again one eigenvalues goes again through zero. The 
squares remain unstable because the other eigen- 
value remains positive. 

Since for 6i ^ the bifurcation point between 
solutal and intermediate squares vanishes, the first 
crossroll pair that emerges from the squares, the 
A1B1A2- and ^ii?ii?2~crossrolls must become dis- 
connected from the square branch. But the oscil- 
latory bifurcation and the other stationary bifur- 
cation remain. For other parameters the scenario 
is simpler. All eigenvalues stay real and only one 
crosses the zero axis. When this is the case, only 
stationary crossrolls appear. 

Note that the pair of eigenvalues responsible 
for the appearance of oscillatory crossrolls and the 
eigenvalue of the stationary crossrolls are not inde- 
pendent of each other. The latter emerges from the 
former. This is a feature that can also be observed 
in the full system of hydrodynamic field equations. 

A bifurcation diagram for slightly different pa- 
rameters and especially larger symmetry-breaking 
terms 6i is shown in Fig. El The stability proper- 
ties are the same as in Fig. ^ They are not dis- 
played here. The bifurcation diagram shows all the 
properties already deducted from the eigenvalues. 
The lower solutal, intermediate, and upper thermal 
branches have fused to a smooth curve that nev- 
ertheless still shows the characteristic upturning at 
the transition from the Soret to the Rayleigh regime. 
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FIG. 8. Bifurcation diagram for Si ^ 0. Symbols iden- 
tify rolls, squares, and stationary crossroUs as in Fig. 
Vertical lines refer to oscillatory crossroUs. The stability 
properties are not displayed here. 

The A1B1A2- and ^ii3ii?2-crossrolls have become 
disconnected from the square branch. The other 
stationary crossroll bifurcation does still exist and 
is preceded by an oscillatory crossroll bifurcation. 
Qualitatively, the picture is exactly as in Fig. ^ 

Before we conclude, we want to take a closer look 
at the dynamics of the oscillatory crossroUs. Fig. IHl 
shows oscillations near the beginning (upper plot) 
and the end (lower plot) of the oscillatory crossroUs 
branch. Comparing to Fig. [21 one sees the same be- 
havior: An oscillation in counterphase with grow- 
ing amplitude, decreasing frequency, and growing 
anharmonicity. 

The details of the entrainment process that lead 
from oscillatory crossroUs to the stationary type 
shall not be discussed here, but they are different 
from what we found in the numerical simulations 
of the full system of field equations jT2] although 
similar transitions have been found in experiments 

We found that the details of the dynamic behav- 
ior of the model system H14|) depends significantly 
on the chosen parameters. Thus, at present we con- 
clude that our model equations ((TH) are capable of 
qualitatively reproducing the bifurcation scenario of 
Fig. ^ and the dynamics of crossroUs in Fig. [2 as 
demonstrated in Figs. |H1 and El 



m 15 - 




t (arb. units) 



FIG. 9. Oscillatory crossroUs in the model for the param- 
eters of Fig. inland = 11.66 (top), fi = 12.14 (bottom). 

5 Conclusion 

We have investigated the Rayleigh-Benard convec- 
tion in binary mixtures with a positive separation 
ratio. We have shown that the observed bifurca- 
tion scenario at small L, involving square patterns, 
rolls, stationary and oscillatory crossroUs arises nat- 
urally out of a simple system of two coupled ampli- 
tude equation systems both consisting of two cubic 
equations describing the dynamic of rolls in x- and 
y-direction respectively. Attempts to extract such 
a model from the basic equations are underway. 
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